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Abstract. Timing analysis of a —95 ks long XMM — Newton observation of the Seyfert 1 galaxy MCG-6-30-15 
is presented. Model-independent tools have been used with the intent of resolving the different components that 
produce the observed flux and spectral variations down to timescales as short as —300 s. We find that the fractional 
variability is possibly due to a variable power law component in the hard band, and a slower thermal component 
at softer energies, but the data do not rule out effects due to the warm absorber. The most relevant result of this 
work is the first detection of a significant enhancement of the fractional variability in the —4.7-5.8 keV energy 
band compared to the underlying continuum. This represents one of the strongest, model-independent, pieces 
of evidence to date for the presence of redshifted relativistic matter accreting into the black hole. During the 
brightest flare of the observation, a soft-to-hard time lag of —600 s is measured. A very significant increase of the 
iron line intensity is observed —3000 s after this flare. In the framework of a disk-corona model and assuming the 
soft-to-hard time lag is due to Comptonization, these findings make it possible to estimate the electron density 
and the dimensions of the flare region and the disk-flare distance. 
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1. Introduction 

According to the commonly accepted paradigm, Seyfert 
galaxies are thought to harbour at their center a super- 
massive black-hole surrounded by a geometrically thin ac- 
cretion disk. The hard power-law component that gener- 
ally dominates the spectral emission above 2 keV is be- 
lieved to arise in a hot plasma corona located above the 
surface of the accretion disk, where optical/UV photons 
from the disk are comptonizcd to X-ray energies. These X- 
rays in turn illuminate the disk, being either "reflected" 
towards the observer or thermalized back into optical/UV 
emission (i.e. the so-called two-phase model, Haardt & 
Maraschi 1991; Haardt & Maraschi 1993; Haardt et al. 
1997). 

Spectroscopic evidence for the two-phase disk-corona 
scenario is seen in X-rays in the forms of a broad, fluo- 
rescent FeK line at 6.4 keV, an FeK edge at <^7 keV, a 
Compton "hump" at >10 keV and a high-energy cutoff at 
E <200 keV (e.g. Nandra & Pounds 1994; Haardt 1997; 
Fabian et al. 2000; Perola et al. 2002). The unprecedented 
large collecting area of XMM — Newton (Jansen et al. 
2001), now offers for the first time the possibility of per- 
forming detailed, spectrally resolved, timing analysis of 
the brightest Seyfert galaxies. These studies potentially 
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offer a method for constraining the emission models and 
their geometry (e.g. Edelson et al. 2002, Poutanen 2001). 

MCG-6-30-15 (z=0.00775) is one of the brightest 
known Seyfert 1 galaxies. It shows a very broad FeK line 
and has remarkable flux and spectral variability proper- 
ties, with doubling times as short as —1000 s. It has been 
observed twice by XMM — Newton: once in June 2000 
for -100 ks (Wilms et al. 2001), and then in July 2001 
for —300 ks (Fabian et al. 2002). In both observations, 
the time-averaged spectral analysis indicated the pres- 
ence of a strong and broad FeK line extending down to 
about 3-4 keV, later detected also by Chandra (Lee et al. 
2002). Detailed model- independent temporal analysis of 
the longer 2001 observation by Fabian et al. (2002) and 
Vaughan et al. (2003a) illustrated well the potential of this 
type of study. The RMS spectrum presented by Fabian et 
al. (2002) showed a decrease in the variability at high ener- 
gies interpreted as due to an interplay between a variable 
power law and a stable reflected component. Vaughan et 
al. (2003a) showed that the power spectral density (PSD) 
of MCG-6-30-15 is very similar to that of Galactic Black 
Hole Candidates (GBHCs), and detected a soft to hard 
time lag of about 200 s. 

Here we present the temporal analysis of the June 
2000 observation, with the data taken from the XMM — 
Newton public archive. The PSD and structure functions 
are given in La Palombara et al. (2002). Those authors 
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show evidence for the absence of any periodicity, and the 
presence of variability down to timescales of ~300 s. We 
perform here further analysis, with particular attention 
to the study of spectral variations and time lags dur- 
ing the whole observation and during a prominent and 
strong flare. The paper is organized as follows. Section 
2 describes the data reduction. In Sect. 3 the results of 
the analysis performed over the whole observation are 
presented, while Sect. 4 presents the results during the 
strongest flare. Thermal Comptonization models predic- 
tions are compared with obtained results in Sect. 5, while 
in Sect. 6 we attempt to give constraints on some physical 
parameters (e.g. density and dimensions of the disk-corona 
system). Alternative models are briefly described in Sect. 
7, followed by conclusions in Sect. 8. 

2. Observations and Data Reduction 

The data were collected during the observation performed 
on June 11 th 2000, with the EPIC MOS1, MOS2 and pn 
cameras (Striider et al. 2001) operated in timing, full- 
frame, and small window modes, respectively. MOS1 tim- 
ing data were excluded from the analysis because it was 
strongly contaminated by the high background. MOS2 
data were also excluded because they suffered from sub- 
stantial pile-up. In the following we thus used only the pn 
data. 

The data were reduced and screened using the SAS 
software v5.3. The analysis of the pn data was performed 
using both single and double events (i.e. pattern <4). An 
a-posteriori check confirmed that the results are identical 
using only single events. The last ~5 ks of the observation 
were excluded from the analysis because strongly contami- 
nated by soft-p + flares. Source counts were extracted from 
a circular region of 45" radius centered on MCG-6-30-15, 
while the background events were taken from two larger 
rectangular regions far from the source. The total use- 
ful exposure time was ^95 ks, yielding a total number of 
—1.2 x 10 6 source counts in the 0.2-10 keV band. The av- 
erage 2-10 keV flux of F 2 -w = 2.4xl0~ n erg cm~ 2 s _1 
is somewhat low for MCG-6-30-15, being only 20% higher 
than during the "deep minimum state" found by Iwasawa 
et al. (1996). Furthermore, during the observation, the 2- 
10 keV flux varied in the range 1.5-4.8x 10 -11 erg cm~ 2 

3. Analysis of the whole observation 

Figure 1 shows the 0.2-2 keV (upper panel), 2-10 keV 
(middle panel) light curves and their ratio (lower panel) as 
a function of time. As already known, the source is highly 
variable with variations up to a factor of —4. The light 
curves show many shots/bursts, with a large prominent 
flare that occurred just before the end of the observation 
at t — 9.2xl0 4 s (see Figure 1, upper panel). Hereinafter 
we will refer to this event as the "Flare" . 

Figure 1 (lower panel) also clearly shows that this 
flux variability is associated to strong spectral changes, 
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Fig. 1. (Upper panel) Soft (0.2-2 keV) and (middle panel) 
hard (2-10 keV) EPIC pn light curves with 100 s time bins. 
(Lower panel) Hardness ratios ( q 2-2fcev ) as a function of time. 
Note that count rates are not corrected for the "live time" 
(~71%) for the pn small window operating mode (Striider et 
al. 2001). 



confirming previous findings on Scyfert galaxies (e.g., 
Mushotzky et al. 1993; Lee et al. 2000; Shih et al. 2002) 
that the spectrum becomes softer as the flux increases. 

3.1. Fractional Variability or "RMS Spectrum" 

To characterize the source variability as a function of en- 
ergy, we used the root mean square variability function, 
also called "RMS spectrum" (for definition and details, 
see the Appendix and Edelson et al. 2002). Similar to 
the more commonly used "excess variance" (Nandra et al. 
1997, Edelson et al. 2001), the RMS spectrum allows the 
quantification of the source fractional variability as a func- 
tion of energy in a model-independent way. The advantage 
of the RMS spectrum (as defined in the Appendix) is that 
its "modified" treatment of the errors allows a finer energy 
binning, thus yielding higher sensitivity to fine spectral 
structures. 

Figure 2 shows the obtained normalized RMS spec- 
trum. It clearly shows three different gross "patterns" : 
a soft, smooth increase of the fractional variability from 
about 0.2 to 0.6 keV; a "plateau" of maximum variability 
between ~0.6-2 keV; a gradual decrease of the variability 
at E >2 keV (of the form oc£T - 16 ). Overall, the gross 
shape and normalization of this RMS spectrum is consis- 
tent with that observed one year later during the 2001 
observation (Fabian ct al. 2002). In particular, the nor- 
malization of the RMS spectrum is slightly less than the 
one from Fabian et al. (2002). This is in agreement with 
the flux-rms correlation found in Vaughan et al (2003a). 

Most importantly, thanks to the finer sampling 
achieved at high energies, we detect for the first time some 
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Fig. 2. RMS spectrum calculated with time bins of 6000 s and 
with energy bins grouped in order to have >350 counts per bin. 
Our choice undersamples the detector energy resolution at low 
energies. Errors are la (see Appendix). 

fine structures above 5 keV, i.e. close to the FeK complex: 
a steep rise of variability in the ^4.7-5.8 keV energy band 
followed by a drop at energies close to 6.4 keV, and a flat- 
tening at energies greater than ^7 keV. To estimate the 
significance of the variability rise in the 4.7-5.8 keV energy 
band, we first fitted the 2-10 keV RMS spectrum with a 
power-law continuum. The fit then improved significantly 
(Ax 2 ~ 15, i.e. >99% confidence in a F-test) when adding 
to the model a gaussian line at 5.3 keV, confirming that 
this excess is highly significant. A detailed discussion and 
interpretation of the source variability is given in Sect. 5.1. 

3.2. Study of the Time Lags 

To search for time lags, we used the cross correlation func- 
tion (CCF; White & Peterson 1994) that we calculated 
using the crosscor package in FTOOLS v5.2 (this com- 
mand implements the DCF, e.g. Edelson & Krolik 1988, 
in the "correl" routine of the "Numerical Recipes" , Press 
et al. 1992). 

The energy intervals for the CCF were chosen so as 
to sample the RMS variability pattern. Light curves were, 
thus, extracted for the following energy bands: Ei (0.2-0.6 
keV), E 2 (0.6-2.2 keV), E 3 (2.2-2.6 keV), E 4 (2.6-4.5 keV), 
E 5 (4.5-6.8 keV), and E 6 (6.8-10 keV). 

Due to the presence of an instrumental switch-off in the 
middle of the observation (see Fig. 1), the CCF study has 
been performed separately, before and after it. Moreover, 
we calculated the CCFs with and without the Flare pe- 
riod. E 2 was used as the energy band of reference for all 
CCFs because the light curve in this band had highest 
statistics. 

Figure 3 shows, as example, the CCFs obtained during 
the 33 ks after the switch-off (which does not include the 
Flare period). The CCFs are fairly wide and show no sig- 
nificant time lag between the different energy bands, with 
an upper limit of about 200 s. The CCF appear also to 
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Fig. 3. Rescaled CCFs calculated in different energy bands 
(Ei:0.2-0.6 keV, E 2 :0.6-2.2 keV, E 3 : 2.2-2.6 keV, E 4 :2.6-4.5 
keV, E 5 :4.5-6.8 keV and E 6 :6.8-10 keV) during the 33 ks fol- 
lowing the instrumental switch-off. Values of CCFs have been 
increased by 1, 0.8, 0.6, 0.4 and 0.2 going from soft to hard 
(top to down) for clarity purposes. 

be somewhat asymmetric, with a marginal excess of cor- 
relation moving from negative to positive time lags as the 
energy bands go from soft to hard (i.e. from top to down 
in Fig. 3). Indicative could also be the possible secondary 
peak in the hardest band corresponding to a time delay 
of ^1.5-2 ks. Similar results have been found, extended to 
higher energies, by the analysis of the simultaneous XMM- 
Newton and RXTE data (J. Wilms, private communica- 
tion). However, in our data the skewness parameter does 
not give any significant deviation from a symmetric profile. 
Better statistics are thus needed for any firmer conclusion 
on this point. CCF analysis of only the first part of the 
observation give the same inconclusive results. 

4. Analysis of the Flare period 

The present observation is characterized by the presence 
of a bright flare, an enlargement of which is shown in Fig. 
4. In studying this event one hopes to isolate the physi- 
cal mechanism(s) responsible for the variability so as to 
constrain emission models and, possibly, generalize these 
findings to other periods when the flares are more elusive. 

4.1. Spectral variability during the Flare 

The light curve hardness ratios shown in Fig. 1 (Sect. 3) 
clearly demonstrate that MCG-6-30-15 exhibits spectral 
variations during the Flare. To quantify this effect, we first 
produced a set of ratios (commonly called "PHA ratios" ) 
between spectra acquired in 8 subsequent time intervals of 
300 s during the Flare (see Figures 4 and 5) and the spec- 
trum accumulated during the whole observation. Figure 
5 shows the results of such an analysis. Before the Flare 
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Fig. 4. Enlargement of the Flare period (see also Fig. 1). The 
light curve is computed between 0.2-10 keV with time bins of 
150 s. Periods during the Flare where PHA ratios have been 
calculated (see Sect. 4.1 and Fig. 5) are labelled from A to H, 
and are 300s long. 



(panel A), the source shows a spectrum close to the aver- 
age one. At the beginning of the Flare, the source shows 
an increase of the soft counts (panel B, C and D) followed 
by a delayed response at higher energies (panel E and F). 
Finally the source gradually returns to its initial spectral 
state (panel G and H). Overall, this result suggests the 
presence of a soft-to-hard time lag. 

To quantify this effect in more details, we performed a 
CCF analysis limited to the Flare period using the same 
energy bands defined in Sect. 3.2. As shown in Fig. 6 (left 
panel), we find that the harder light curves systematically 
lag the softer one by - 80 s, 260 s, 310 s, 340 s, and 480 
s, respectively. This soft-to-hard lag is even clearer from 
Fig. 6 (right panel) where the peak time lag is reported as 
a function of energy 1 . The time lag clearly increases with 
energy (n ag <~ 0.071 x E; where E is in eV and the time 
lag is in s). 

This result is consistent with the time lag (~ 200 s) 
that was found during the 2001 long (-300 ks) XMM - 
Newton observation by Vaughan et al. (2003a). In that 
case, however, the time lag was found only on timescales 
longer than ~10 4 s. This was interpreted as evidence for a 
drop of coherence on shorter timescales. The fact that we 
detect it in a single Flare event (during which the source 
is likely to be dominated by a single active region, thus 
with maximum coherence) , coupled with the null result on 
the whole observation (Sect. 3.2) confirms and strengthens 
this hypothesis. 



1 To calculate the errors associated to the time lag, we fol- 
lowed two different methods: CCF peaks were first "x 2 "-fitted 
with a Gaussian function and the errors were then estimated 
by assuming i) a Ax 2 ~ 4, and ii) -^=, where a is the FWHM 
of the Gaussian and N is the number of time bins considered 
in the fit. Most conservative values obtained following either 
method were used and reported in Fig. 6. 



4.2. FeK reverberation in response to the Flare? 

To study the FcK line variations during the Flare period 
in a model-independent way, we extracted the source spec- 
trum in three time intervals (1000 s long) before, during 
and after the Flare (we will refer to these time intervals as 
the PreFlare, Flare and PostFlare periods). The PHA ra- 
tios between the PostFlare/PreFlare periods highlights a 
strong line increase that took place during the PostFlare 
period. During the Flare only a softening of the contin- 
uum is visible. To further investigate the variations, we 
performed a spectral analysis on a set of eight subsequent 
time intervals of ~1000 s each around the Flare. We fit- 
ted these data using a simple power-law model (r ~1.75- 
1.95) plus a Gaussian line at ~ 6.4 keV. Figure 7 (Right 
Panel) clearly shows that only upper-limits on the EW (of 
~ 100-200 eV) were obtained for the FeK line before and 
during the Flare. On the other hand, the line recorded af- 
ter the Flare (periods 6 and 7 in Fig. 7) was remarkably 
strong (EW~680±200 eV and 280±150 eV) and broad (a 
~0.39±0.20 and 0.19±0.12 keV, respectively). It should be 
noted that a similar effect was found for the absolute flux 
of the FcK line, with A FcK < 6.5 x 10" 5 ; = l.S 



-0.3 



x 10" 



< 1.6 x 10~ 5 photons cm" 2 s~ 2 at periods 5, 6 and 7 in 
Fig. 7 (right), respectively. More complex models for the 
line (e.g. a DISKLINE model; Fabian et al 1989) have not 
been applied given the low signal-to-noise ratio of these 
1000 s spectra. 

It is interesting to note that such spectral changes, de- 
tected here, are qualitatively similar to those found by 
ASCA in 1994 during another remarkably bright flare 
(Negoro et al. 1999; Negoro et al. 2000) 

5. Implications of the observed spectral variability 

The X-ray emission mechanism operating in Seyfert galax- 
ies is usually thought to be inverse Compton on seed soft 
photons that are iteratively up-scattered by a hot elec- 
tron corona (e.g. Sunyaev & Titarchuk 1980; Haardt & 
Maraschi 1991; Haardt, Maraschi & Ghisellini 1997). We 
test this hypothesis against the spectral variations mea- 
sured during the XMM — Newton observation. 

5.1. Constraints from the long (> 5 ks) timescale 
variability: Continuum and FeK line 

As noted earlier (Sect. 3.1), the RMS spectrum allows, 
for the first time, the detection of fine structures at the 
energy of the FeK line complex. In particular, the drop 
in variability at ~6.4 keV is consistent with a fraction 
of the FeKa line being produced in outer (less variable) 
parts of the accretion disk (see also Fabian et al. 2002; 
Lee et al. 2002). The variability excess located between 
4.7 and 5.8 keV instead implies the presence of an intense, 
broad spectral component that varies more than the un- 
derlying continuum. This excess could be due to a broad, 
redshifted iron line component produced by fluorescence 
from an accretion disk around either a Schwarzschild or 
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Fig. 5. MCG-6-30-15 spectral variations during the Flare. PHA ratios are calculated by dividing spectra acquired in subsequent 
time intervals of 300 s of about 4500 counts (see Fig. 4) and the average spectrum of the entire observation. This figure illustrates 
the gradual softening, followed by a hardening and subsequent decrease of the spectrum during the Flare. 




Fig. 6. (Left panel) CCFs between light curves extracted during the Flare. Energy bands are as in Fig. 3, and values have been 
increased by 2, 1.6, 1.2, 0.8 and 0.4 going from soft/top to hard/down for clarity purposes. The dashed vertical line indicates 
a zero lag value. For each energy band, the differences between the peaks of CCF distributions and the dashed vertical line 
indicate the time lags. (Right panel) Time lag (in s) vs. energy (in eV) is shown. Peak time lag values and associated (68% 
confidence) errors have been estimated through a x 2 -fitting method (see footnote n. 1). 



a Kerr black hole. If so, the emitting iron should have a 
bulk velocity of about 0.3-0.4 c. Another, maybe more ex- 
treme, possibility could be that the excess of variability 
is associated to absorption, rather than emission, from al- 
most fully ionized iron falling into the black hole (along 
our line of sight) at about ^0.2-0.4 c (see e.g. Nandra et al. 
1999, Ruszkowski & Fabian 2000, Longinotti et al. 2003). 
Whatever its precise origin, the excess of fractional vari- 
ability detected here between ~4.7 and 5.8 keV represents 
strong (because model-independent) evidence for the pres- 
ence of relativistic matter accreting (because redshifted) 
onto a black hole. 

In the soft band, the sharp drop of variability at E 
<0.6 keV is likely due to the rise of a sharp soft-excess 



component clearly seen in the source-average spectrum 
(e.g. Fabian et al. 2002; Pounds & Reeves 2002). This fea- 
ture, best fitted by a black-body component of average 
temperature kT ~100 eV, is consistent with either intrin- 
sic or reprocessed thermal emission from the less variable 
optically thick accretion disk. On the reprocessed emission 
hypothesis, the data do not exclude the presence of soft 
X-ray disklines produced from the inner, more variable, 
disk regions which may contribute between ~ 0.6-1.0 keV 
(Branduardi-Raymont et al. 2001; Lee et al. 2001; Sako 
et al. 2003; Turner et al. 2003). Another contribution to 
the variability in this band could come from the warm ab- 
sorber that has been observed to vary on timescales down 
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Fig. 7. (Left panel) Flare/PreFlare and PostFlare/PreFlare PHA ratios. (Right panel) FeKa line equivalent width (in eV, left 
y axis) calculated in 1000 s spectra, over-plotted on light curve (right y axis) from Fig. 4. Errors are at 90% confidence level for 
one interesting parameter. 



to 10 4 s (Otani et al. 1996), and therefore should be ob- 
servable by the RMS spectrum. 



In the hard band, the drop of variability could be due 
to different intcrplaying phenomena. Phenomcnologically, 
it could be explained as a single power-law spectrum that 
varies with a pivot point at E >10 keV, as the analysis of 
the hardness ratio indicates (the spectral index softens as 
the flux increases). Physically, such behavior is consistent 
with what is expected if the coronal optical depth varies 
while the dissipated luminosity remains constant, and it 
is just what is expected for the disk corona models as long 
as the spectral index rr 2 _io] is Zz>2 (Haardt et al. 1997). 
Another possibility could be to have a two-component con- 
tinuum model like the one proposed by Shih et al. (2002), 
Fabian & Vaughan (2003) and Miniutti et al. (2003), 
which foresees one steep power law rigidly variable plus 
a flatter and almost constant reflection-dominated com- 
ponent. In this model, the reflection component is pro- 
duced very near the black hole and is heavily affected by 
light bending. As a result, the emission line can be strong 
down to ^3-4 keV (Fabian et al. 2002), and could reduce 
the continuum variability in a complex pattern. Another 
possible explanation for the complex FeK line versus con- 
tinuum variability is the ionized reflection model proposed 
by Ballantyne et al. (2003) which requires two distinct 
reflectors: one from the inner accretion disk to build up 
the red wing and the other from the outer disk to fit the 
blue peak. Alternatively, some contribution to the vari- 
ability decrease at hard energies could also come from the 
Comptonization process itself. In fact, higher energy pho- 
tons, after a great number of interactions (i.e. at larger 
distances and in different physical states), may display a 
smoother variability than do softer photons. 



5.2. Constraints from the short (~ 100 s) timescales 
variability: Continuum and FeK line 

The spectral variability during the strong (AF~2 in 
~1000 s) Flare is characterized by a soft-to-hard time lag 
of a few hundred seconds (see Fig. 6, right panel). If the 
Flare time lag, Ti ag ~600 s, is interpreted as the delay re- 
quired by Comptonization to upscatter soft photons up 
to 10 keV (see also §6.1), then it is comparable with, or 
even larger than, the minimum time scale (~300 s) of the 
variability. 

The remarkably strong (EW ~680 eV) FeK line is de- 
tected with a delay of ^3000 s after the Flare (Fig. 7, 
right panel), which may be indicative of a disk reverber- 
ation in response to the Flare. This indicates that simu- 
lations should also incorporate this geometric delay when 
computing the expected spectral variations. Moreover, we 
note that, just after the Flare, the FeK complex is present 
up to ~7.0 keV. This suggests that the Flare may have 
been responsible for ionizing the disk and/or illuminating 
its approaching (blue-shifted) parts. 



6. Mapping the physical parameters of the 
disk-corona region 

In the framework of a hot corona model scenario, and 
assuming that the lag observed during the Flare is due 
to Comptonization, it is possible to derive a number of 
physical and geometrical parameters of the disk-corona 
system, as follows. 

6.1. Dimensions and density of the Flare region 

Assuming that the electron temperature in the emission 
region, best constrained by the high-energy BeppoSAX 
data (Guainazzi et al. 1999; Fabian et al. 2002), is in 
the range ~100-300 keV, seed soft photons emerging 
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from the disk will need about 5-10 interactions 2 (A^ nt ) 
to be up-scattered from ^200 cV to 10 keV. Thus, the 
mean free path of 10 keV photons (Am) is linked to the 
Comptonization timescale (tcompt) through the equation: 
Am = tcompt x c = T '^ t c ~ 2-4x 10 12 cm, where r lag 
is the time lag (^600 s) between the 0.2 and 10 keV emis- 
sion and c is the speed of light. Assuming a random walk 
motion and N ln t=10, this implies a displacement from 
their original production region of <~ VMnt Aio ~ 6 x 10 12 
cm. Therefore it is reasonable to assume that, when the 
bulk of 10 keV photons are produced, the Flare dimensions 
(^Fiaro) should be of the order of <~ 6 x 10 12 cm. For val- 
ues of the optical depth (t) ranging between r <~ 0.1 — 1, 
a typical range for Seyfert galaxies including MCG-6- 
30-15 (e.g. Petrucci & Dadina 2003), this yields also a 
rough estimate of the Flare region electron density (p e -) : 
Pe~ = a- T R Fl — ~ ^ 25 x 10 10 cm -3 , where ut is the 
Thomson cross-section. 

The dimensions of the Flare region (i?Fiarc ~6x 10 12 
cm) estimated above correspond to ~ 40 gravitational 
radii (r g ) assuming the mass of the black hole in MCG-6- 
30-15 to be Mbh ~ 1O 6 M (e.g. Vaughan et al. 2003a). 
This value is consistent with the limit implied by causality 
arguments for the factor 2 variation in ^1000 s during the 
Flare (Fig. 4), which gives i?Fiaro^ c x t ~ 3 x 10 13 cm ~ 
200 r g . 

6.2. Measure of the disk-corona distance 

The spectral analysis performed during and after the Flare 
(Fig. 7) clearly demonstrates a significant increase of the 
FeK line intensity ~3000 s after the Flare peak. Assuming 
that the line is produced in the accretion disk and in re- 
sponse to the Flare, it is possible to estimate the distance 
between the Flare emission region and the disk. To per- 
form this accurately, one should know when precisely the 
intensity of the line started to increase, the transfer func- 
tion of the disk and the physical condition of the Flare 
region. These are all unknowns but, if ever significant, 
they would introduce delays. Thus the estimate given be- 
low shall be considered only an upper limit. 

Since the hard radiation ionizing the iron on the disk 
surface lags the Flare peak by about ~600 s (§4.1, Fig. 6), 
the upper-limit on the distance between the Flare region 
and the disk (A; F _ d ) is given by: A/ F _ d < 3000 2 ' 600 x c = 
4 x 10 13 cm <~ 240 r g , also consistent with the marginal 
presence of a second peak in the CCFs at higher energies 
discussed in §3.2, Fig. 3. 



2 We estimated the number of interactions that a photon 
undergoes to be up-scattered from ~200 eV to 10 keV by iter- 
ating the equation hv ~ ■y^hv that links the average photon 
energy before (hv) and after (hv ) an inverse-Compton scat- 
tering event. In this equation, 7 e is the electron Lorentz factor 
given by the equation: E e — m e c 2 (7 e — 1), and we assumed that 
the electron energy (E e ) is equal to the high-energy cut-off of 
the X-ray spectrum. 



7. Other models for the observed spectral 
variability 

Instead of considering models with a single corona (as 
in the earliest two-phase models) one could also consider 
patchy corona models (e.g. Haardt et al. 1994), i.e. mod- 
els in which several active regions are present, with (more 
or less) independent physical (i.e. spectral) evolution. In 
this scenario, the short timescale variability is expected 
to be dominated by the evolution of the individual active 
regions. On long timescales, these models must be able to 
reproduce the observed broad-band spectra that clearly 
suggest thermal Comptonization. However, the single flare 
event could be dominated by different, non-thermal, mech- 
anisms. In fact, there are problems in maintaining a ther- 
mal plasma in rapidly varying sources, since it may not 
have time to thermalize (Ghisellini et al. 1993). 

Alternative mechanisms could for example involve 
magnetic reconnection above the disk surface (Haardt et 
al. 1994; Stern et al. 1995; Poutanen & Svensson 1996; 
Bcloborodov 1999a, b; di Matteo et al. 1999). Such mecha- 
nisms can introduce changes in the energy dissipation rate 
and/or in the geometry of the flaring regions (e.g. their 
distance from the disk) and could therefore reproduce the 
soft-to-hard spectral evolution found in the present paper 
(e.g. Poutanen & Fabian 1999). Electron-positron pairs 
could also contribute significantly to the coronal optical 
depth during the strongest flare periods. Indeed, using the 
standard equation for the compactness I = ^ (assum- 
ing L = L 2 -io keV during the Flare, and R = 3 x 10 13 ), 
we estimated a value of I ~10 during the Flare period 
(corresponding to a pair optical depth ~0.1, Stern et al. 
1995), indicating that, at least during the most compact 
flares, the pair contribution may not be entirely negligible 
for the source spectral evolution (Haardt et al. 1997). 

If such alternative mechanisms are indeed respon- 
sible for the time lag observed during the Flare, 
the Comptonization timescale must be correspondingly 
shorter. In this scenario, the estimate of the electron den- 
sity of the Flare region given in §6.1 should be considered 
a lower limit. 

Finally, we anticipate that in the future more pre- 
cise timing analysis will be able to discriminate among 
these alternative scenarios. In fact, models in which (short- 
and long-term) spectral variations are driven by only 
Comptonization and/or disk-corona reverberation effects 
should exhibit only soft-to-hard time lags, contrary to 
other models that may have zero lag and/or hard-to-soft 
time lags. 

8. Conclusions 

The model-independent timing analysis of a <~ 95 ks long 
XMM - Newton observation of MCG-6-30-15 presented 
in this work has shown that: 

(1) The fractional variability around the FeK line en- 
ergy is "resolved" with unprecedented detail into two dif- 
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fcrcnt components: a narrow component of lower variabil- 
ity at E ~6.4 keV consistent with its being produced by 
the outer parts of the accretion disk and a broad com- 
ponent in the 4.7-5.8 keV energy band, significantly more 
variable than the continuum itself (Ballantyne et al. 2003). 
Whether this second component is due to emission or ab- 
sorption is unclear, but it is likely associated to a highly 
variable and redshifted iron line that probes the inner re- 
gions of the accretion disk. 

(2) The reduced variability of the source continuum 
at soft (E <600 eV) energies is consistent with a gradual 
onset of soft thermal emission from the outer parts of the 
accretion disk. The reduced variability at higher (E >2 
keV) energies is possibly due to a power law having a 
pivot point at E >10 keV, or to a gradually increasing 
contribution from a reflection component. 

(3) During a remarkably strong Flare, a significant 
soft-to-hard time lag of ^600 s is detected. It is the first 
time such a time lag has been detected in MCG-6-30- 
15 on such a short timescale. We find it to be consis- 
tent with a delay due to Comptonization up-scattering. 
On this hypothesis, we are able to infer estimates of 
the electron density and dimensions of the Flare region: 
p e - ~ 2.5 - 25 x 10 10 cur 3 and i? Flarc - 6 x 10 12 cm 
<~ 40 r g , respectively. 

(4) Ratios of spectra taken during and after the Flare 
show evidence that the FeK line intensity has varied very 
significantly, and may have done so in response to the 
Flare with a ^1500-3000 s lag. In the simplest scenario, 
in which the Flare is produced in the corona and the FeK 
line is reflected from the disk, the observed lag allows us 
to estimate the disk-corona distance as £4 x 10 13 cm ~ 
240 r g . 

(5) The present short-term-variability results well il- 
lustrate the potential offered by a high-throughput instru- 
ment like XMM — Newton for obtaining constraints on 
theoretical emission models, and map the inner regions of 
AGNs. 
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Appendix A: Definition of the RMS spectrum 

We used an innovative fractional variability function, called "RMS spectrum" . It allows a reduction of the associated 
errors, and therefore a finer energy binning, when a sufficiently high number of counts is available. If on this condition 
of high statistics per bin, one can consider negligible the Poisson noise, the standard fractional variability amplitude 
(F V ar, e.g. Nandra et al. 1997, Edelson et al. 2002), can be rewritten as: 



RMS(E) 



N 



[x(E,AU)- < x(E) >p 
N — 1 



< x(E) > 



(A.l) 



where x(E,Ati) is the source count rate in the energy band E during the time interval At,, < x{E) > is the mean 
source count rate at energy E during the whole observation, N is the number of time bins (i.e., N= ExposureTime y rp-^ 
RMS spectrum characterizes the source fractional variability as a function of the energy. Compared to the more widely 
used F var , the RMS(E) has the advantage that its associated errors can be calculated analytically by the standard 
propagation equation, yielding 



CRMS 



where a. 



N r 



dRMS 
dx(E, AU 



■dx(E,Ati) 



N 



J2[(x(E,AU)- < x{E) >)<r x(EAU) f 
\ »=i 



< x{E) > 



N 



(A.2) 



i=l 



x(E.Ati) is the error associated to x(E,Ati), and where we assumed that the error associated to < x(E) > is 
negligible compared to (J x (e ,au)- 

This assessment of the error is less "conservative" than using the more widely used fractional variability equation 
(Fig. A.l, left panel) since it requires that the energy and/or the time bins used in the calculation be chosen so as 
to have a negligible random Poisson noise. Under these conditions, the error given in Eq. A.2 (Fig. A.l, right panel) 
is also very similar to the one calculated recently by Vaughan et al. (2003b) using Monte Carlo simulations (see Fig. 
A.l, middle panel and see Eq. B.3 in Vaughan et al. 2003b). Its greatest advantage is however that in cases like 
the present one of large signal-to-noise, it does reduce significantly the errors and thus allows a much finer sampling 
of the energy and/or time scales involved (cf. Edelson et al. 2002 for the overestimation of the errors caused by a 
conservative approach in the calculation of the standard fractional variability function) . Furthermore, since the errors 
and the RMS data points are formaly independent (unlike those calculated by Vaughan et al. 2003b) we could estimate 
the significance of the variability structures fitting directly the RMS spectrum (see §3.1). In particular, in the present 
data, we chose energy and time bins in order to have at least 350 counts per bin, so as to exclude any Poisson noise 
(which we found to be relevant only for <200 counts per bin). Finally, we checked that the results obtained using the 



1000 2000 
energy (eV) 



1000 2000 
energy (eV) 



1000 2000 
energy (eV) 



Fig. A.l. RMS spectra calculated with the equations used in Edelson et al. 2002 (Left panel), in Vaughan et al. 2003b (Middle 
panel) and in this work (Right panel). The main difference of the RMS functions is their error estimates, as follows: Edelson 
et al. (2002) used a conservative approach in the calculation of the standard fractional variability function, leading to an 
overestimation of the errors; Vaughan et al. (2003b) estimate the errors with Monte-Carlo simulations; the formula presented 
here is a good approximation of the Vaughan et al. results if Poisson noise is negligible (i.e. if > 300 counts/bin). The advantage 
of the equation presented by Vaughan et al. is that it is always valid, also in the case of few counts/bin, but its disadvantage is 
that it cannot be formally be fitted to estimate the significance of structures, since its errors are not independent. The advantage 
of the RMS equation presented in this work is that the errors are independent, so fine structures can be fitted. The disadvantage 
of this method is that it is applicable only if the statistics are good enough to ignore Poisson noise (> 300 counts/bins). 
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RMS spectrum are entirely consistent with those obtained using a standard fractional variability formula. But thanks 
to the lower error bars, we could detect the fine structures with increased significance (see text). 
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